Baryons and baryonic matter in the large Nc and heavy quark limits 
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This paper explores properties of baryons and finite density baryonic matter in an artificial world 
in which Nc, the number of colors, is large and the quarks of all species are degenerate and much 
larger than Aqcd- It has long been known that in large Nc QCD, baryons composed entirely 
of heavy quarks are accurately described in the mean-field approximation. However, the detailed 
properties of baryons in the combined large Nc and heavy quark limits have not been fully explored. 
Here some basic properties of baryons are computed using a variational approach. At leading order 
in both the large Nc and heavy quark expansions the baryon mass is shown to be Mbaryon ~ 
NcMq (l — 0.05426 Qg) where fis = NcCtg. The baryon form factor is also computed. Baryonic 
matter, the analog of nuclear matter in this artificial world, should also be well described in the mean- 
field approximation. In the special case where all baryons have an identical spin fiavor structure, 
it is shown that in the formal heavy quark and large Nc limit interactions between baryons are 
strictly repulsive at low densities. The energy per baryon is computed in this limit and found to be 
exponentially small. It is shown that when the restriction to baryons with an identical spin-flavor 
structure is dropped, a phase of baryonic matter exists with a density of 2A^/ times that for the 
restricted case but with the same energy (where Nf is the number of degenerate flavors). It is shown 
that this phase is at least metastable. 

PACS numbers: 11.15.Pg,12.39.Hg,14.20.-c 



I. INTRODUCTION 

One of the central theoretical problems in nuclear 
physics is the understanding of the properties of cold nu- 
clear matter directly from QCD. However, this problem is 
quite difficult since, unlike the study of QCD at non-zero 
temperatures, QCD at non-zero baryon chemical poten- 
tial cannot be studied using standard lattice techniques 
due to the presence of a fermion sign problemfT. In- 
deed, from the point of view of the functional integral 
of QCD it is not even understood why the system at 
zero temperature remains in the vacuum state for non- 
zero chemical potentials below the critical value for pro- 
ducing nuclear matter; this is the so-called Silver Blaze 
problem[2]. Given this situation it would be useful to try 
to gain some intuition into the problem via the study of 
limiting cases of QCD for which the baryon matter prob- 
lem (the generalization of the nuclear matter problem) 
becomes tractable. 

One obvious approach for seeking tractable limits, is to 
start with the large Nc limit of QCD. It has long been rec- 
ognized that certain aspects of QCD greatly simplify in 
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the large Nc limit [3j. At the fundamental level of quarks 
and gluons, the large Nc theory is dominated by planar 
graphs, while at the phenomenological level, the zero- 
baryon sector of the theory becomes a weakly coupled 
theory of mesons and glueballs. Baryons are more com- 
plicated. One way to proceed is to consider a theory in 
which one flavor of quark is in the two-color antisym- 
metric representation and other flavors are in the color 
fundamental representations ,4, 5j. Such a theory has the 
virtue that 3-quark baryons exist at any Nc- However, 
such an approach badly breaks flavor symmetry and thus 
is an unnatural generalization of the Nc = S world[5]. 
An alternative approach, pioneered by WittenlT], is to 
respect flavor symmetry and consider baryons which con- 
tain Nc quarks and deal with the combinatoric complica- 
tions that arise due to this. We will follow the 't Hooft- 
Witten large Nc limit here in which gluons are in the 
adjoint representation and quarks are in the fundamen- 
tal. 

Unfortunately, despite its simpliflcations large Nc 
QCD is not tractable in general. Thus, it cannot be used 
directly to gain insight about baryonic matter. To pro- 
ceed further, additional simpliflcations are needed. One 
approach is to restrict attention to large Nc QCD in 1+1 
dimensions — the 't Hooft model [5. Meson spectroscopy 
in this model is tractable: it can be obtained from so- 
lutions to a simple integral equation which can be eval- 
uated numerically. Underlying this simpliflcation is the 
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fact that the ghion degrees of freedom are non-dynamical 
in 1+1 degrees; there are electric components of the gluon 
field strength but by construction no magnetic ones and 
thus giuons cannot propagate. Recently, there have been 
significant advances in the understanding of baryon mat- 
ter in the 't Hooft model [9llT0|. 

One might similarly hope to find tractable regimes 
in 3-1-1 dimensional large QCD where the dynam- 
ics associated with the gluon sector is absent or highly 
restricted. This occurs if the characteristic momentum 
flowing through giuons is much larger than Aqcd . In 
this case, the running coupling is small and dynamical 
gluon effects are suppressed. One way to achieve this is in 
the regime of very high density or equivalently chemical 
potential, /i . In that case, /i provides a large scale which 
can push the typical momentum flowing through giuons 
well above Aqcd ■ Large Nc analyses of the high regime 
have been done both for the usual 't Hooft large Nc limit 
(with quarks in the fundamental representation) [llj and 
for QCD^g — a variant in which quarks are in the two in- 
dex anti-symmetric representation — which represents a 
different extrapolation from Nc = 2> [H]. For the case of 
the 't Hooft large Nc limit it was shown that a DGR[T3] 
instability occurs, indicating that the ground state can- 
not be of the BCS type but rather breaks translational 
symmetry [TT]. However, in QCD^g the DGR instability 
can be shown not to occur and it is plausible that the 
ground state is of the BCS typed^- The differences be- 
tween these two examples illustrate a critical point. In 
going to extreme limits to gain tractability, one may well 
push the system to a regime which differs in important 
qualitative ways from QCD at Nc = 3 in less extreme 
conditions. The existence of a DGR instability in the 't 
Hooft large Nc limit might be viewed as an artifact of the 
large Nc limit. The regime of interest involves a double 
limit: fi ^ oo, Nc — >■ oo. It is not obvious that these two 
limits commute and, in fact, they do not: as stressed in 
ref. jl4] the behavior depends sensitively on whether the 
the large Nc limit or the high density limit is taken first. 

This paper focuses on another regime of large Nc QCD 
for where dynamical gluon effects are suppressed: the 
case where the quark masses (which are taken to be de- 
generate) are much larger than Aqcd- In this case the 
characteristic momentum flowing through dynamical giu- 
ons is typically a^Mg. The quarks, being very heavy, 
move nonrelativistically in this regime. The giuons be- 
come nondynamical: they give rise to a color-Coulomb 
interaction between the quarks, and effects of propagat- 
ing giuons are suppressed. It is thus not surprising that 
this regime provides a tractable setting for the study of 
both single baryons and baryonic matter. 

The combined heavy quark and large Nc limit has a 
long history. Indeed, Witten's classic paper on the prop- 
erties of baryons in large Nc QCD is based on an analy- 
sis of baryons in this combined limit; it is first argued 
that this limit is accurately described in a mean-field 
approximation From this, the Nc scaling of various 
baryon properties were deduced. Finally, Witten argued 



that these properties were generic and remained valid 
even away from the heavy quark limit. Witten's reason- 
ing shows the power of using tractable regimes to learn 
rather deep things about large Nc QCD more generally. 

This paper explores baryons and cold baryonic mat- 
ter in the combined large Nc and heavy quark limits. 
The problem of baryons is in a general sense quite well 
understood and the logic underlying Witten's mean-field 
approach to baryons is extremely well known. However, 
while the formulation is over three decades old, to the 
best of our knowledge the equations have never been 
solved numerically. It is important to do so if for no other 
reason than completeness — it remains one of the only 
regimes of QCD which is solvable. There is also some 
phenomenological utility in doing this — the Vli, baryon 
(composed of three bottom quarks) might be regarded 
as being relatively near this in this limit. While the J7f, 
has not been observed to date, there is little doubt about 
its existence in the real world P^. 

The second thrust of this paper is the problem of cold 
baryonic matter. We treat the problem of infinite baryon 
matter using the mean-field approach suggested by Wit- 
ten. The general solution to this problem has yet to 
be solved and it is not immediately clear whether bary- 
onic matter in the combined limit saturates (i.e., is self- 
bound). However, we can show that there exists a phase 
at low density which consists of weakly interacting nu- 
cleons. In this regime the nucleus- nucleus potential is 
due to the Pauli principle and strictly repulsive; how- 
ever, the repulsion is exponentially small. We can com- 
pute the equation of state for cold matter in this phase. 
It is unknown whether this phase is absolutely stable. It 
is known, however, to be at least metastable with a very 
long lifetime. The question of absolute stability will need 
to be investigated in future work. 



II. THE BARYON IN THE COMBINED LARGE 
Nc AND HEAVY QUARK LIMITS 

A. The mean-field energy functional 

The general large Nc scaling rules for baryons were 
originally deduced by Witten [7]. The analysis in ref. [7] 
begins by considering the special case where the quarks 
constituting the baryon are heavy. The formalism for the 
case of heavy quarks is rigorous. Somewhat more heuris- 
tic arguments were then invoked to suggest that the Nc 
scaling results for the case of heavy quarks should also 
hold when the quarks are light. This allows for the de- 
termination of the large Nc scaling rules for baryon in 
the phenomenologically interesting regime of nearly chi- 
ral quarks. Here we will stick to the heavy quark regime. 
The analysis of this regime in ref. [7] was based largely on 
physical reasoning. In this subsection a somewhat more 
formal derivation will be given. Aspects of this more for- 
mal treatment help clarify the nature of the mean-field 
approximation for the baryon. Moreover, it will prove 
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extremely useful in formulating the baryonic matter cal- 
culation in the following section. 

The starting point of the analysis in ref. [7] is that in 
the heavy quark regime, quark-antiquark pairs are sup- 
pressed and the quarks behave nonrelativistically. Color 
magnetic couplings are also suppressed since the color 
magnetic moment scales as the inverse of the quark mass. 
Moreover, gluon-gluon interaction effects are suppressed 
since the characteristic size scale of the system is small 
and thus ag at this scale is small (albeit only logarith- 
mically). Thus for heavy quarks, a baryon is accurately 
described as Nc quarks in a color singlet state interact- 
ing with each other through a color Coulomb interaction. 
The nonrelativistic Hamiltonian for this system in second 
quantized notation is 



\ C \C 

O^s ij i' j' 



2Mc 



Nc 2 2 
with cts = OisNc , 
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where = |^ is the strong coupling constant (at the 
appropriate scale for this problem), otg (the ' t Hooft 
coupling) is independent of Nc^ Mq is the (common) 
quark mass, q are nonrelativistic fermion field operators 
with spin and flavor indices suppressed (and implicitly 
summed over), ^ are generators of SU(Afc) and i, j, i', j' 
are color indices; there is an implict summation over c 
and 

The Hilbert space of physical baryon states consists 
of color singlet states with Nc quarks. The color-singlet 
nature of the state imposes anti-symmetry with regard 
to color with each color represented once and only once. 
This fact implies that matrix elements of the Hamiltonian 
in Eq. ([T]) only receive non-vanishing contributions for 
terms in which the color of each annihilation operator q 
is identical to the color of one of the creation operators, 
q\ in that term. Thus in the color Coulomb term there 
are two types of terms which could contribute: those with 
i = j, i' = i' and those with i ~ j' , i' ^ j : 



^3^ ^3^> 4i^4'i'r')m'{r')qiir) 



Nc 2 2 

^ /■ 3 3 q\ {r)ql (r')ft' (r)% (r') 

Nc 2 2 ' ^ 



(2) 



The first term's contribution is at most order Nc- Note 
that it can only contribute for those generators with 
nonzero diagonals. There are only A^^ — 1 of these (each 
with two nonzero diagonal matrix elements) and thus the 
sum over c contributes a factor of at most Nc which is 
canceled by the \/Nc factor out front. Since the lead- 
ing order energy is order Nc this term may be dropped. 
Moreover, one can use the fact ^u'^iH = 2{1—Sii' /Nc) 



to simplify the color Coulomb term to be 

The anti-symmetry in color implies that the many- 
body states must be symmetric remaining degrees of 
freedom — space, spin and flavor. Thus for n-quark states 
in the physical state one can suppress the color degrees 
of freedom and treat the quarks as (spin-1/2) bosons. 
Introduce bosonic operators, Qfs{r), which satisfy 

[Qfs{r),Qf,s'{P)]^0 

[Ql{r),Q\,Ar')]=0 (4) 

[Qfs{r1,Q\,Ar')] = 6fj'5,,s'SHr- Sr') . 

Since the color factors in Eq. ([3| reduce to simple sum- 
ming over all of the quarks, one can rewrite the nonrela- 
tivistic Hamiltonian as 



H= I d\Q\r) [Mq- 
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(5) 



where we have dropped terms suppressed in 1 /Nc and for 
simplicity have suppressed explicit reference to the spin 
and flavor. In what follows in this section, we will use 
the bosonic version of the theory. 

Note that at leading order in a heavy quark expansion, 
as we have here, neither the spin nor flavor of the quarks 
is dynamically relevant -the only role they play is to give 
the degeneracy of the baryon state. Accordingly, in our 
initial treatment it will be assumed that all quarks are 
in a single spin and flavor state and the spin and flavor 
degrees of freedom will be ignored throughout. The issue 
of degeneracies due to spin and flavor degrees of freedom 
will be discussed at the end of this section. In general, 
even with the neglect of the spin and flavor degrees of 
freedom, this is a very complicated quantum mechani- 
cal many-body problem. However, the problem greatly 
simplifies at large Nc- 

The physics underlying this is quite simple. The in- 
teractions between any two quarks makes a small contri- 
bution to the energy — of order l/Nc- However, provided 
that the pairs contribute coherently as in a mean-field 
state, the fact that there are of order Nc pairs, leads 
to a contribution of order Nc and, thus, contributes at 
leading order. Thus the key issue is the coherence of the 
state. It is apparent from the symmetry under exchange 
of particles that the optimal coherent state will have each 
particle in the same quantum state. 
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where ■(/'(f) is the optimal spatial wave function. The 
expectation value of the energy in states of this form is 
given by 



1 



_ * f-UL q 

2 1 1^' _ ^'1 



(7) 



with J = 1 



where we have dropped subleading terms in l/Nc- This 
is the form used in the mean- field analysis of ref . [7 ; our 
expression differs slightly in convention using the stan- 
dard normalization for g in which ag = g^/(47r). It is 
convenient to rewrite N^ag as N^as which makes clear 
that all three terms on the right-hand side of Eq.([7| are 
proportional to Nc . 

Minimizing (iplHlip) with respect to tplr) gives the 
leading order result for the mass of the baryon. We 
denote this state IV')!!!!!! and the spatial wave function 
"0111111(^0- The wave function ip(f) is obtained by mini- 
mizing the energy. Starting with Eq. ([?]), a straightfor- 
ward application of the calculus of variations yields an 
equation for ipmini'r)- 

^ = V'miiilrO 

(8) 

where e is a Lagrange multiplier enforcing the normaliza- 
tion. 



B. The nature of the mean-field approximation 

States of the form of Eq. ^ define a variational sub- 
space. Clearly the nature of variational calculations is 
such that we can use it to provide a bound on the en- 
ergy. However, a question arises as to how well does this 
state actually approximate the true ground state of the 
Hamiltonian in Eq. ([5| at large Nc- As will be discussed 
in this subsection, generically one expects 10)111111 to be 
a poor approximation to the true ground state, |V')Gxact, 
even in the large Nc limit: the overlap between |'0)cxact 
and 1 0;) mill does not approach unity as Nc — >■ 00. Instead 
one has 



1 



0{N"c) 



(9) 



P « 1. 



In practice, it is quite possible that | exact (010)1 
However, despite the fact that |0)niin is a poor approx- 
imation to the true quantum state, min ('01-^^10) mm is 
a good approximation to the baryon mass Mbaryon = 
exact (01^^10) exact: fractiou errors are of order 1/Nc and 
thus go to zero at large Nc- This may seem paradoxical 
at first sight but it has a very simple physical origin: in 
the exact ground state at large Nc, the probability that 



all of the quarks have the wave function ^min is small 
while the probability is very high that most (that is, all 
except a fraction of order l/Nc ) of the quarks have this 
wave function and thus dominate the energy functional. 

To understand why this is so, it useful to decompose 
the Hamiltonian into a one-body mean-field Hamiltonian 
and correction terms: 



H = H„,i + AH 
AH=V - + C 



(10) 



with 
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The form of the one-body mean-field Hamiltonian 
is seen in the one-body Schrodinger equation of 
Eq. It is easy to see that min(0|V^|0')min = 

min(0|14if|0)min ^ and C is choscu so that 

mill (01 AT? I t/i) mill = 0. The eigenstates of the one-body 
mean- field Hamiltonian, 77mf, provide a complete basis 
of states. By construction, |0)min is the ground state 
of -ffmf- Excited states of -ffmf are n-particle-n-holc ex- 
citations relative to |0)min- In the large Nc limit, the 
excitation energy of a generic n-particle-n-hole state is 
independent of Nc- 

While 10) mill is the ground state of 7?mf it is clearly 
not the ground state of H . The presence of AH in H 
introduces correlations in the form of n-particle-n-hole 
state contributions to the ground state. At first blush it 
may seem reasonable to regard AH as a perturbation and 
use perturbation theory to compute these correlations. In 
implementing such a calculation it is useful to understand 
the Nc scaling of matrix elements of AH between various 
n-particle-n-hole states (including the 0-particle-O-hole 



state 



i): 



((n ± l)p - (n ± l)h| Ai?|np - nh) - Nc 
((n ± 2)p - (n ± 2)h| Ai?|np - nh) - TV," 
((n ± m)p — (n ± m)h| Ai? |np 



nh) =0 for TO > 2 

(12) 



At first order in perturbation theory, there is no correc- 
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tion due to AH. At second order, 



AM, 



bar yon 



nelp-lh " " 

^ |(n|Ai?|^)„,i„|2 



(13) 



Tie2p-2/i " 



where the superscript nif indicates that these are the en- 
ergies in the unperturbed mean-field Hamiltonian. From 
Eq. (12 1 it is apparent that the contribution from the 



one-particle-one-hole states are of order 1/Nc while the 
contributions from the two-particle-two-hole states are of 
order N^. Since the mean- field calculation gives a mass 
of order , this correction is down by a factor of order 

iv-i 

The correction to the energy in perturbation theory 
is parametrically small compared to the leading order 
contribution. This is what one would expect if the sys- 
tem was in the regime of validity of perturbation theory. 
However, if one looks at the correction to the state, one 
sees that the dominant (2-p-2-h) contribution is given by 



E I 



lAi/r 
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mf 



mf 



(14) 



Clearly, low-order perturbation theory cannot be justified 
on parametric grounds at large Nc] corrections to the 
state are generically not small. 

The failure of low-order perturbation theory in AH 
seen in Eq. ( 14 1 implies that the perturbative compu- 
tation of the mass shift in Eq. ( 13 ) is unreliable. This 



presents a major challenge in any attempt to compute 
subleading correction to the mass in the 1 /Nc expansion; 
these include contributions due to perturbations due to 
AH. However, for the purposes of extracting the leading 
order behavior, it is not necessary to compute the cor- 
rection accurately: the relevant issue is the Nc scaling of 
the contribution to the mass due to the existence of AH. 

Fortunately, perturbation theory does capture the 
1 /Nc scaling properly. To see why, one can imagine doing 
the following artificial calculation. Consider the eigen- 
spectrum of 



Hx = Hn,i + XAH 



(15) 



where A is a parameter. The actual problem of interest 
corresponds to A = 1. In this artificial problem AAfbaryon 
is a function of both A and 1/Nc. Let us consider what 
happens for small A. In this case we can justify pertur- 
bation theory in XAH: 



AAfba 



CXD 

E 

3=0 



(16) 



Using the scaling rules in Eq. ( 12 ), it is easy to show that 
each of the coefficients rrij scales as at leading order 



in the \/Nc expansion. Expanding the coefRcients nij in 
a series of 1/Nc yields a double expansion: 



AAU 



baryon 



EE^^'^. 

j=0 fe=0 



rrijk 



(17) 



where the coefficients rrijk are independent of A and Nc- 
Providing that the double expansion is well behaved, one 
can switch the order of summation yielding 



AMbaryon ^^Nc^Mj{X) with 

OO 

k=0 



(18) 



In effect, in introducing Alj we have resummed the A ex- 
pansion to get the contribution to all orders. The leading 
order term in the 1/iVc expansion thus scales as Nc for 
any value of A. Extrapolating to A = 1 shows that cor- 
rections to Mbaryon duc to the prcseucc of AH are generi- 
cally of absolute order N^ which is of relative order 1/Nc 
compared to the leading order contribution. The argu- 
ment given here does have a loop-hole in that it depends 
on the technical assumption that the double expansion 
is well behaved. We have no reason to doubt that this is 
the case. 

Thus, the mean-field approach advocated by WittenjT] 
appears to be valid as a technique for obtaining the lead- 
ing order value of the mass in a 1/Nc expansion. How- 
ever, it does not give the leading order wave function. 
Moreover, the 1/Nc expansion is not equivalent to a per- 
turbative expansion in AH . 

The apparent discrepancy between mean-field theory 
providing an accurate computation of the energy while 
yielding an inaccurate wave function is easy to under- 
stand. At a physical level it is apparent in the exact 
ground state at large N^ the probability that all of the 
quarks have the wave function ?Ainin is small; this explains 
the breakdown of mean-field theory for the wave func- 
tion. At the same time only a small fraction of the quarks 
(of order 1 /Nc ) are not in the mean-field single-particle 
level. Thus the mean-field quarks dominate the energy 
functional. From this physical picture, it is immediately 
apparent that one can accurately compute observables 
beyond the energy in mean-field theory. In particular, 
any n-body operator with n <^ Nc can be calculated; the 
expected error will be of relative order n/Nc. 

One final comment about the nature of this mean-field 
approximation. The analysis above was in the framework 
of the effective bosonic theory. Of course, the result is 
the same if one works in the original fermionic theory. In 
that case, instead of using the coherent state of Eq. ^ 
one would use a Slater determinant. From the form of the 
interaction term in Eq. (|3|, it is apparent that the term 
which contributes is the exchange interaction and not 
the direct interaction. This means that the contribution 
comes from the Fock term (i.e., the exchange interaction) 
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and not the Hartree term {i.e., the direct interaction). 
Thus, the description in ref. [7] of the approximation as 
a Hartree approximation is something of a misnomer. It 
is important to keep this in mind when generahzing the 
treatment to the problem of baryonic matter, as is done 
later in this paper. 



C. Parametric dependence 

Before explicitly implementing a variational calcula- 
tion to find the baryon mass, it is useful to understand the 
parametric dependence of the mass on Nc, as and Mq. 
The Nc dependence is trivial: the mass is clearly propor- 
tional to Nc. More generally, the mass of the baryon can 
be seen to be given by Mbaryon = NcMq (l — const a^). 
This can be shown via a simple scaling argument. Note 
that the right-hand side of Eq.([7]) has three terms: a 
quark mass term, a quark kinetic term, and an attrac- 
tive interaction term due to color Coulomb interactions. 
One can ask how these three terms scale with the spa- 
tial size of the wave function? Introduce a length scale 
parameter R with 



(19) 



where foiv) is a reference normalized functional form 
which characterizes the sliapc of the wave function; y 
is dimensionless and the size of the wave function is thus 
characterized by R. 

{i^\H\,p)^Nc{MQ+t + v) 



M I AT , Tjfo) ' ^ifo) 

v,/o(y)P 



r(/o) ^ / d'y 



V(/o)= / d^yd'y 



(20) 
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t and V are the kinetic and interaction terms, and T and 
V are dimensionless functionals of the reference function 
/q. Suppose ipo is the wave function that minimizes the 
energy. In that case we know that dii{iJj\H\ip) = 0. From 
this we deduce that 



R = 



t 



1 2r(/o) 

asM V(/o) 

2 ■ 



4r(/o) 

li^lV^) = A4aryo„ - iV,A/Q (l 
.;;2 V\fo) 



(21) 



Mr(/o) 



Thus NcMQCtj ^j-Yf„) can be regarded as the binding en- 
ergy of the heavy quarks. Note that an overall dilation 
of /o is innocuous: it changes the equilibrium value of 
R, T and V while keeping {i/j\H\iIj) unchanged. A direct 



comparison of Eq. ([8| with Eqs. (20) and (21 ) shows that 



t + 2v 



4 r 



(22) 



Thus the problem of finding the mass of the baryons 
in this limit reduced to determining a single function /o 
from which T and V and, hence, the mass of the baryon 
can be computed. The functional form for /o does not 
depend on cts or Mq and thus neither do T and V. 



D. A variational calculation 

Equation ([7| is a nonlinear integro-differential equa- 
tion. One standard strategy in obtaining an approximate 
numerical solution to such equations is iterative. An al- 
ternative approach which is pursued here is to consider 
a variational space of trial wave functions. A wave func- 
tion in such a space is specified by some finite set of pa- 
rameters 6i , 62 . . . , 6n • Evaluating the energy in Eq. ([7| , 
using states in the variational space and minimizing with 
respect to the parameters 61 , 62 ■ ■ ■ , ^'n will yield an opti- 
mal wave function for the space of trial wave functions. 
If the space is rich enough to include the exact solution of 
Eq. ([8|, then the wave function obtained by minimizing 
with respect to 61 , 62 ■ • ■ , bn will necessarily be a solu- 
tion of Eq. ([s]). If the trial space does not contain the 
exact solution but does contain states with a very large 
overlap with the exact solution, one expects that solu- 
tion obtained via minimization in the trial space to ap- 
proximate the exact solution with high accuracy. The 
virtue of this approach is that one finds explicit analytic 
expressions for the wave functions — albeit only approx- 
imate ones. Moreover, there is considerable experience 
indicating that even quite modest variational spaces can 
yield remarkably accurate results for energies. 

One can represent a function in the trial class of 
normalized dimensionless functions as /o''^'^^'' "''"' (y) so 
that the wave function is given by ipi^^^'''^'---^^) [r) — 
R-^/^f^''''^^'-''"\f/R). The expectation value of the en- 
ergy in this state is given by 



EiR, 61, 62, . . . bn) = {4,('>iM,-K)\H\^p«^,.b. 



,.fc„)\ ^ 



NcImq 



R^Mr 



(61,62, ...6„)n 



R 



(23) 



The optimal values for the coefficients are obtained from 
minimizing E{R, &i, &21 • • ■ bn) with respect to R and the 
hi. Extremizing i?(i?, 5i, 62, . . . &„) yields the following 
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set of coupled equations: 



562 



(24) 



aiogrgr-"--^")) 

56„ 



= 2 



91ogV(/, 



db„ 



For any particular class of functions, Eq. (24) is simply 
a set of algebraic equations and is amenable to standard 
numerical solutions. It is easy to verify that a particular 
solution is a (local) minimum as opposed to a maximum 
or saddle point. 

To implement this scheme in practice one needs to 
choose a space of functions which is rich enough to closely 
approximate the exact solution. There is a strong theo- 
retical prejudice that the exact solution should be spher- 
ically symmetric. While it is certainly possible for mean- 
field solutions to break symmetries, in the present cir- 
cumstance it seems clear that a spherical shape both min- 
imizes the kinetic energy and maximizes the magnitude 
of attractive color Coulomb potential energy. Thus, our 
trial class of functions consists of spherically symmetric 
functions. 

One useful way to construct a trial class of functions 
is to begin with a complete set of spherically symmetric 
functions, {(/)o(?;), 0i(y), 02(2/), ■ • •}, which satisfy an or- 
thonormality condition, dy 4>*{—)4>l.{y) = 5j i^. The 
exact form of /o can always be written as a normalized 
superposition of the . A truncation of the sum in the 
superposition yields a useful trial class of functions: 



Ei 



(t>o{v) 



(25) 

This form is viable provided that the shape of an exact 
solution to Eq. Q has a nonvanishing overlap with 0o. 
By truncating the expansion at a different value of n one 
can test numerical convergence. A particularly useful 
complete set of functions for this construction is 



2/V(2j + l)! \^ 



1 d 



2dy) 



ye 



(26) 

these are the s-wave eigenfunctions of the harmonic oscil- 
lator (up to an irrelevant scale factor) and a sign conven- 
tion with 0j(O) positive. These functions have the virtue 
that an analytic expression for energy in terms of the bj 
can be computed straightforwardly. Moreover, with this 
trial class of functions the energy converges quite rapidly 
as additional states are included into the sum. 

We have minimized the energy function using trun- 
cations with up to seven basis states, (00, ■ • ■ ^e)- The 
results for T and V are given in Table I. It is apparent 



TABLE I: This table shows the numerical convergence of the 
coefficients T and V defined in Eq. ( 20 1 for the minimum en- 



ergy configuration computed in a truncated basis. The func 
tions are given in Eq. ( 26 1 
the binding energy of the heavy quarks 



The combination |^ determines 



basis 
functions 


4r 


r 


V 


00 


0.053052 


1.5000 


0.56419 


00, 01, 02 


0.054198 


1.4293 


0.55666 


00,01, ... ,04 


0.054252 


1.4067 


0.55250 


00,01, • • • ,06 


0.054256 


1.3941 


0.55006 



that the quantity ^ which determines the binding en- 
ergy of the heavy quarks is well converged and is given by 
0.05426 up to four significant figures. Thus in the heavy 
quark and large Nc limits the mass of the baryon is well 
approximated by 



« N,Mq (1 



0.05426 5^1 



(27) 



The wave function at this level of truncation is deter- 
mined by the coefficients 5i . . . fcg which, to good numer- 
ical approximation, are given by 



5i « -0.0476828 , 62 
64 ~ 0.0136048 , 65 ~ 



0.0815798, 63 R 
-0.00159968, fog 



-0.00786028 , 
0.00283365 

(28) 



with the scale factor R given by 

1 2r(/o) 5.06907 



R 



asM V(/o) 



(29) 



Combing these yields the following approximate wave 
function: 

V'mi„(r) « 0.0876207(a^MQ)3/2 exp (-0.0389173f2) 
X (0.755925 - 0.00856005f2 _^ o.000289408f4 

1— I o cocnoin— 8-8 



-1.6598710 
with f = rd^M, 



10-10 



3.6419310- 



-13-12N 



(30) 



The wave function is plotted in Fig. [TJ 

It is clear from Table I that the variational procedure 
converges rapidly yielding an accurate expression for the 
mass. Of course, it is the nature of variational approxi- 
mations that energies are computed more accurately than 
the wave functions: first order difference between the ap- 
proximate and exact wave functions do not contribute to 
the energy. A useful way to test the accuracy of the varia- 
tional wave function itself is to compare the left-hand side 
of Eq. (30) with -0; for the exact wave function the dif- 
ference is zero. This comparison is made in Fig. [T] where 
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FIG. 1: 

against r in the solid curve. V is given in units of {osiy-iq 



The variational wave function of Eq. (|30|) is plotted 

3/2 

1 



and r is given in units of {asMq)~^ . The short-dashed curve 
represents the left-hand side of Eq. ([8| using the variational 
wave function of Eq. (30 1. For the exact wave function, the 
two coincide. 



the short-dashed curve represents the left-hand side of 
Eq. ( 30 ) using the variational wave function and the solid 

Note that overall 
The largest 



curve is the variational wave function, 
the two curves are quite close everywhere, 
difference between the two is at the origin and amounts 
to a 1% discrepancy. Elsewhere the discrepancy is quite 
small for the entire region where the wave function is 
substantial. For typical observables this is the region of 
relevance and thus we conclude that we can accurately 
compute typical matrix elements with high precision — at 
least for a world with and Mq very large. 

Note that the form of all of the wave functions in the 
variational space used here drop off at very large dis- 
tances quite rapidly — like a power law times a Gaussian. 
As will be shown below, the exact wave function drops 
off as a power law times a decaying exponential, which 
is much slower than the states in the variational space. 
This means that at very long distances the variational 
wave function is much smaller than the actual solution 
of Eq. ([t]). Fortunately the variational space used is rich 
enough so that, although the wave function is inaccurate 
at asymptotically large distances, it is quite accurate over 
the domain containing the vast bulk of the baryon's prob- 
ability. Hence the variational wave function is capable of 
describing virtually all of the standard baryon observ- 
ables of interest such as the mass, the baryon charge ra- 
dius or the form factor in the range where ~ (Mga)^. 
However, if one focuses on observables which are domi- 
nated by the behavior at a very large distance, the wave 
function in Eq. (30) is inadequate. 



Fortunately, it is possible to determine the long dis- 
tance behavior of the wave function accurately. One 
can start with the variational wave function Eq. ( 30 ) 



and determine the form of the correct asymptotic of the 
solution Eq. ([?]). One can match this onto the varia- 
tional wave function in a domain where both the asymp- 
totic form and the variational wave function are accu- 
rate. The asymptotic form is found easily: for large r. 



IV-min(r')|^ 
\r—r' I 



tical to the Schrodinger equation for a Coulomb poten- 
tial. Of course the solutions are not hydrogenic — even 
asymptotically — since the eigenenergy, e, does not corre- 
spond to those of hydrogen atom. Using this asymptotic 
form of the potential and solving the Schrodinger equa- 
tion yields an asymptotic wave function of 

V'Asyl'^) = (Mq 5,) ^Asy(r) with 



r 



U is the confluent hypergeometric function of the sec- 
ond kind. The constant, /3, can be fixed by matching 
this asymptotic form to the variational solution. The 
matching needs to be done in a regime where r is large 
enough so that the corrections to the asymptotic form 
of Eq. (0 are negligibly small while simultaneously be- 
ing sufficiently small so that errors in the variational 
wavefuntion due to its wrong asymptotic behavior are 
also negligible. By comparing the r dependence of both 
the variational and asymptotic wave functions we found 
that these conditions were well satisfied in the region 
7 < f < 11. We chose to match at f = 9 and found 



f3 w .118 



(32) 



We estimate the error in the determination of /3 to be 
less than .5%. 



E. The baryon form factor 

One of the most important quantities characterizing a 
baryon is its form factor. For a nonrelativistic system 
such as the one being considered here it is defined by 

G{q') ^ {p'\pi6)\P) 
2 



(33) 



N, 



and Eq. (7) asymptotically is iden- 



where \p) is a baryon state of good momentum defined 
with the normalization {p'\p) = (2Tr)^ S^{p — p') and q is 
the quark field operator. G is essentially the relativistic 
Sachs form factor G^; in the regime where the momen- 
tum transfer is much less than the baryon mass. For 
the equivalent bosonic theory that we are using here, the 
density operator may be written as Pb{x) = Q (^Q(^) ^ 

One typically expects that in mean-field models that 
break translational invariance, the form factor is given by 
the Fourier transform of density in the center of mass: 



(34) 
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It is useful to verify that Eq. ( 34 ) does indeed follow from 



the mean- field approximation, and that it is justified at 
large Nc for heavy quarks. In doing so one can determine 
the parametric dependence of the corrections to this re- 
sult. 

As noted in Subsec. |II B[ one expects the mean-field 
formalism to be valid at large Nc for few-body opera- 
tors. As p is a one-body operator it appears that mean- 
field theory should be appropriate. There is a subtlety, 
however, in that the mean-field state breaks translational 
invariance and thus is not an eigenstate of momentum; 
accordingly it cannot be used directly in Eq. (34). In or- 



der to compute the form factor it is necessary to decom- 
pose the mean-field state as an integral over momentum 
eigenstates 



I V') mill ^ 
A{p)\p)n,f : 

P : 



(27r)3' 



(35) 



(f'xQ^x) (-«v) Q{x) 



|p)inf is the (normalized) state of good momenta obtained 

by projecting from the mean-field state and P is the mo- 
mentum operator. It is straightforward to determine the 
parametric dependence of A(j)) on Mq and fisi 



A{p) = 



A 



( \ 



(36) 



where the functional form of A is independent of Mq 
and Ois- 

Consider J d?'xe^^'^^ia{'il>\pb[x)tp)mia- On the one 
hand, using Eq. (6]) and the definition of this is easily 
seen to be given by 



dPxe'^-^min{i^\pb{x)'^)min= / d^X e*«-^|V'mi„(^) | 



(37) 

On the other hand, using Eq. (35) and the fact that 
p{x) — e^^'^ p{0)e~^^''^ , it can be written as 



(fxe'"^'^^in{ll^\pb{x)\'ip)nun = 

|V e^i^^+i'-i^-'A*{f)A{p) ^i{p\MO)\P).ni 



(27r)3 (27r)3 
2^ I d^P 



(27r)3 



G.mCz') / Si: A 



^A*{p)A{p + q) 
(2^)3- (l)i(p + 9(v^«sA/Q) 



(38) 



where Gmf is the baryon form factor defined in Eq. ( 33 ) 




0.0 0.2 0.4 0.6 0. 



1.0 1.2 1.4 



FIG. 2: The baryon form factor in the large Nc and heavy 
quark limits. In the plot, is given in units of {a.sMqf' 



. Note that 



dimensionless form introduced in Eq. ( |36 

at large Nc, with q independent of Nc, <f(\ 
goes to 0; dropping it in the last form of Eq. (38) yields 
Gmiiq^) since the integral of \A\'^ is unity. Thus, at 
large Nc we have Gn,f{q^) = / d^a; e'«' Vin(V'l/56(^)|V')min 
which immediately yields Eq. (34). 



At finite Nc there are correction terms to Eq. ( 34 ) since 
(^^/NcCtsMQ^ 0. It is easy to see that parametri- 

cally they are formally of relative order q^ (TVcSsMq) ^ . 
The interesting kinematic region is for ~ Q^s-^q since 
this is the characteristic scale over which the form factor 
varies. Thus, in this region there is a correction of rela- 
tive order 1 /Nc to Eq. ( 34 ) from the approximation given 
above. There are other corrections of the same order due 
to the fact that the mean-field wave function is not exact. 



We have evaluated the form factor using Eq. ( 34 ) and 
the wave function in Eq. (30). The result is plotted in 
Fig. |2] It is of interest to consider derivatives of the form 
factor with respect to q^ evaluated at = since these 
are related to moments of the baryon density distribu- 
tion: 



) baryon 6 



) baryon — 60 



dG{q^) 
dq^ 

CPG{q^) 



d{q^y 



21.5 
890. 



(39) 



where the numerical values were obtained by the mean- 



field form factor using the wave function in Eq. ( 30 ) and 



evaluated using the mean- field state \'4>)m\-a and A is the 



thus are accurate at leading order in the double expansion 
in l/Nc and Aqcd/Mq. 



F. Degeneracy 

Due to the spin and fiavor degrees of freedom there are 
a large number of distinct species of baryons and, in this 
limit, the masses of many of these degenerate. Recall that 
baryon mass is insensitive to either the spin or the fiavor 
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of the quarks making up the baryon. The reason for this 
is clear — by hypothesis the mass of the quarks are degen- 
erate and in the heavy quark hmit the spin of the quarks 
decouple from the other degrees of freedom in the prob- 
lem. Moreover as discussed above in the mean-field state 
the spatial wave functions are all identical s- waves. Com- 
bining this with the anti-symmetry from Fermi statistics 
and the anti-symmetry associated with the color degree 
of freedom implies that, at leading order in the double 
expansion, the degeneracy of the ground state baryon 
is given by the number of distinct symmetric spin and 
isospin configurations that can be constructed from Nc 
quarks. 

There are 27V/ possible spin-isospin states for each 
quark. The symmetry under exchange implies that each 
state in the degenerate multiplet is fully specified by the 
number of quarks in the baryon, with each of the 2Nf 
possibilities subject to the constraint that the total num- 
ber of quarks is Nc- Thus the degeneracy is given by the 
following bionomial coefficient 



(40) 



Of course, this degeneracy is split. Firstly if the quark 
masses are not precisely degenerate there is a splitting 
due to the mass differences. However, even in the exact 
fiavor symmetric limit there is a splitting due to color hy- 
perfine splittings. Characteristically the maximum split- 
ting of the full multiplet is characteristically given by 
AMmax C(sNcXqcd/Mq. Notc that in this formalism 
the heavy quark limit is taken prior to the large Nc limit 
and thus the total splitting remains small in the com- 
bined limit. 



G. The value of as 

The calculations detailed above treat as as a constant. 
Of course, in reality there is a scale anomaly. Thus 
the coupling constant runs, albeit logarithmically. Using 
standard renormalization group analysis one sees that at 
one loop and at large Nc 




I2n 



llLog (^) 



(41) 



where Q is the momentum scale of the process and Aqcd 
is the QCD scale. Thus, in the formalism used above, 
is to be interpreted as being evaluated at the momentum 
scale of the problem. 

Of course, the question of what constitutes the opti- 
mal choice of scale is not completely well posed. The 
effective value of the coupling differs for different parts 
of the wave function. However, over the region of the 
wave function which makes the dominant contributions, 
as changes by only a small amount since the running is 



logarithmic and the inverse length scale of the wave func- 
tion is much smaller than KqCD. Thus, the ambiguity 
is of relatively small importance provided that the SsMg 
is high enough. The inclusion of higher order corrections 
in d-s would further reduce the ambiguity. 

The relevant momentum scale of this problem is dsMg, 
which is the inverse of the characteristic length scale of 
the wave function. This yields a transcendental equation 
for OLs'. 



Utt 



llLog (52 



's A 



2 

QCD 



(42) 



The analysis in this paper is based on the regime in 



which 



Mq 



> 1 and as < 1- 



These conditions 



might be satisfied to reasonably good approximation even 
well away from the extreme large mass limit in which 

Log ( .2 ^ >• 1; in those cases it is necessary to solve 

VlflCD / 

Eq. ( 42 ) numerically in order to find as , given the value 
of Mq . Using this value of fig in Eq. (271 yields the 
baryon mass. 

/ \ 

However, in the regime where Log ( . ^ ^ 1 , it 

V '^QCD / 

is possible to estimate analytically as- The key point is 
that in this regime. Log (Sg) /Log ( ) ^ 1; thus we 

can rewrite Eq. ( 42 ) as 

127r 



■^QCD 



as 
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(Log( 



Mi 



127r 



llLog 





^^QCD 



+ Log ((52) 



Log (5^) 



^^QCD / 



(43) 



where we have dropped terms of relative order 



Log^ ((52) /Log^ ^ ^2'^ j . Using the leading order result 
to estimate the value (5^ in the subleading term yields 
127r 



llLog ( 



^QCD 



127r 



(Log(^: 



Log Log 



Mr 



(44) 



Log2( 



Mi 



"-QCD 



It is apparent that the second term in Eq. ( 44 1 is indeed 



suppressed compared to the leading one in the regime 
Log I X2~^ ^1. In this regime one can use Eq. (44) to 

V '^QCD / 1 

find (5s which, in combination with Eq. (27), yields the 
baryon mass. 



III. BARYONIC MATTER IN THE COMBINED 
LARGE AND HEAVY QUARK LIMITS 

Baryonic matter is the analog of infinite nuclear mat- 
ter in this artificial world with many colors and heavy 
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quarks. In describing baryonic matter in this regime it is 
necessary to develop a formalism which correctly encodes 
the leading orders in both the large and heavy quark 
expansions. In ref. [7j Witten assumed implicitly that 
the argument justifying mean-field theory for the single- 
nucleon case also holds for the many-nucleon problem. In 
particular, ref. [7j assumed that the appropriate mean- 
field wave function for a two-baryon scattering-state is 
a Slater determinant composed of two distinct clusters, 
each composed of Nc quarks in identical spatial wave 
functions coupled to a color singlet. Before proceed- 
ing, it is probably worth remarking that although this 
is quite plausible, it is not trivially obvious. The multi- 
baryon problem is more complicated than the analogous 
single-baryon case since the bosonic substitution made in 
Subsec. II A no longer holds. However, the qualitative 
reasoning underlying this assumption is quite plausible 
and we will assume it to be correct for the remainder 
of this paper. In particular, it will be assumed that, as 
with the single-nucleon case, such mean-field calculations 
will accurately give the mass (up to corrections of rela- 
tive order l/N^) while not accurately describing the wave 
fimction. 

In the analysis of the single-baryon problem, the spin- 
fiavor structure of the baryon was irrelevant — provided 
that the state was symmetric. The heavy-quark limit 
leads to a decoupling of the spin degree of freedom while 
the degeneracy of the quark masses leads to a fiavor 
symmetry. Thus the only possible effect of the spin- 
fiavor symmetry is due to the effect of the Pauli prin- 
ciple. In the single-baryon case, anti-symmetry of the 
wave function with respect to color means that the space- 
spin-fiavor part of the wave function is symmetric. Since 
the mean-field treatment necessarily leads to a symmet- 
ric spatial wave function, the spin-flavor state is auto- 
matically symmetric. However, in a certain strict sense, 
mean-field theory of the Hartree-Fock type applies only 
to a subset of such baryons: namely states in which all of 
the quarks are an identical spin-fiavor state. These could 
be states where each quark is in well-defined projections 
in spin and flavor space separately or one in which spin 
and fiavor are correlated as, for example, in a hedgehog. 
In other spin-flavor conflgurations the state cannot be 
described as a single Slater determinant since there are 
correlations between spin, flavor and color. Thus, one 
should envision the myriad of states in the single-baryon 
system as a two-part problem: first, one uses a standard 
Hartree-Fock analysis based on a single Slater determi- 
nant to compute the energy. One then subsequently can 
make transformations on top of the single Slater deter- 
minant state which do not alter the energy (for heavy 
quarks) to consider states with spin-flavor correlations. 

In this section, we study the energy density of cold 
baryonic matter using mean-field techniques, which as 
noted above, are assumed to hold in the combined large 
Nc and heavy quark limits. Following Witten, we will 
assume that the state can be written as a single Slater 
determinant and that any effects of spin-fiavor correla- 



tions can by obtained by acting on top of a single Slater 
determinant without changing its energy. In general, 
even with this assumption, the problem is rather com- 
plicated. However, it greatly simplifies at low density. 
We show that there exists a phase which is repulsive and 
has an interaction energy per baryon which is exponen- 
tially small at low densities. This phase is shown to be at 
least metastable. This phase can be the true ground state 
depending on whether baryonic matter in this limit satu- 
rates. While we leave this question to future research, we 
will briefly consider the effects of 1/iVj, corrections and 
show that a phase in which baryonic matter saturates 
is likely to emerge at next-to-leading order in 1/Nc but 
that density and binding energy per baryon of saturated 
matter both go to zero as Nc goes to infinity. 



A. A simplified problem 

As noted above, for the single-baryon case the only way 
to construct a single Slater determinant while imposing 
a color singlet structure is to have all the quarks in an 
identical spin-flavor state. The problem under consid- 
eration in this subsection is to determine what happens 
when there is a non-zero density of identical baryons in 
this class and no other baryons in the system. To state 
the problem more precisely, we seek the ground state of 
a system with a fixed density of quarks each with the 
same spin-fiavor state. We will turn to the more general 
problem without this restriction later. It is worth not- 
ing that the present problem is not well posed in QCD 
with finite quark masses since the spin quantum number 
is not conserved. However, in the limit where Mq — ;> oo, 
the spin degree of freedom decouples and spin becomes 
an emergent symmetry with an associated conservation 
law; in this limit the problem is well posed. 

By fixing each quark to have an identical spin-fiavor 
state — even quarks associated with different baryons — 
the problem simplifies considerably: the spatial wave 
function for two quarks in the systems with identical color 
needs to be anti-symmetric. Moreover, at the mean-field 
level, by construction, there are two and only two distinct 
wave functions for the two quarks of the same color; anti- 
symmetry then implies that they are orthogonal. Thus 
for a system with a total baryon number, A^bary, greater 
than unity we can consider a color singlet state of the 
following form: 



n n / d'''V',(r>L(^') 



j=l a=l 



(45) 



where a indicates color, the subscript sf indicates a par- 
ticular spin-fiavor state and j distinguishes between the 
various single particle states; the wave-functions ipj are 
taken to be orthonormal: / d'^rip*{r)ipk{r) = 6jk- To 
make the infinite matter problem concrete, one should 
envision keeping the number of baryons, A'bary: large and 
finite and then imposing boundary conditions restricting 
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all of the wave functions to have support only within 
some volume, V. The limit, — > cx), iVbary — >• oo with 
Nhaiy/V held fixed is then taken at the end of the prob- 



lem. The many-body state, |$), can then be inserted to 
into the nonrelativistic reduction of the QCD Hamilto- 
nian yielding the following energy functional: 



2^cAs E E (^^^ -1) / 

j=l k=l 



2Mq 



\r' - r'l 



(46) 



subject to the constraint / d^rtpj{r)'fpk{r) — Sjk- The 
first term represents the dominant energy arising from 
the mass of the heavy quarks. The second term may be 
interpreted as the interaction energy of quarks within one 
baryon summed over all baryons and the third term as 
the interaction energy between baryons. Note, however, 
in general the second term differs from the interaction 
energy between quarks in an isolated baryon since the 
quark wave functions rearrange themselves due to the 
presence of the additional baryons. 

The next step is to find the minimum energy config- 
uration subject to the constraints. It is generally be- 
lieved that in the large Nc limit baryonic matter will 
form a crystal. The underlying reason for this is clear 
if one thinks about baryonic matter using a simple pic- 
ture of baryons interacting via potentials: baryons are 
heavy and their interactions strong. In such a regime 
one generically expects baryons to find the minimum of 
the potential created by the others and sit there. The 
kinetic energy of the baryon — which is what typically 
fights against crystallization — scales as 1/Nc and is thus 
suppressed. At a more fundamental level the large Nc 
limit is, in an important sense, classical: the expecta- 
tion value of a product of color singlet operators is the 
product of the expectation values of the operators. Thus, 
quantum correlations of the physical degrees of freedom 
are suppressed. However, spatial correlations are not; ex- 
pectation values of operators can vary in space. To mini- 
mize energy, the system will generally exploit its freedom 
to build in spatial variations yielding either crystals or 
amorphous solids. For a wide class of systems, crystals 
have lower energy than amorphous solids. The general 
behavior of baryonic matter forming a crystal is seen ex- 
plicitly in Skyrme-type models |16j which are designed to 
reproduce the large scaling behavior of QCD. More- 
over, the heavy quark limit also leads to heavy baryons 
and should not affect the qualitative result that crys- 
talline phase is expected. Indeed the heavy quark limit 
should act to reinforce this expectation. 

To compute the energy density, we will need an ex- 
plicit ansatz for the crystalline form. If a low density 
phase of baryonic matter exists, the appropriate ansatz 



which minimizes the energy is easy to find. As seen from 
Sec. [njthe wave function for a single baryon drops off ex- 
ponentially at long distance. Moreover Eq. ( 46 1 implies 
that the interaction energy between baryons depends on 
the product of the two wave functions at the same spa- 
tial point. If the wave functions of interacting baryons 
in a low density system are qualitatively similar to free 
baryons, one expects that the interaction energy between 
two baryons in the system to drop exponentially with the 
distance between baryons. Effects of the Pauli principle 
will also be exponentially small. Thus, in a regular crys- 
talline configuration at low densities all interactions are 
expected to be exponentially small. Moreover, nearest- 
neighbor interactions although exponentially small are 
also exponentially larger than next-to-nearest-neighbor 
interactions. There are two possibilities to consider: the 
case where the nearest neighbor interactions are attrac- 
tive and the case where they are repulsive. For the at- 
tractive case, the low density matter is not stable and the 
system will collapse; low density baryonic matter does 
not exist. However, if the interaction is repulsive — as we 
shall show below it is — then low density baryonic mat- 
ter exists at least as a metastable phase provided that 
there is an external pressure to prevent the system from 
expanding. Moreover, the crystal structure which mini- 
mizes the energy density for fixed baryon density is the 
one which maximizes the nearest-neighbor distance at 
fixed density. 

Let us denote the nearest neighbor distance as d. By 
dimensional analysis, it is related to the baryon density 
by p = Jj- where c is a constant that depends on the 
crystal structure. It is conventional in crystallography to 
describe structures in terms of their atomic packing frac- 
tion P which gives the maximum fraction of the volume 
filled by rigid non-overlapping spheres of fixed radius ar- 
ranged in the given crystalline configuration. It is related 
to the constant c by c = 4^. Thus, 



1/3 



(47) 



and the configuration with the largest nearest-neighbor 
spacing for fixed baryon density is the one with the 
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largest value for the atomic packing-factor, P. It has 
been known since Gauss that the maximum value for P 
in any crystalline configuration is P = which occurs 

for both hexagonal close-packed (HCP) lattices and for 
face-centered cubic (FCC) lattices. Thus, provided that 
a stable low density phase exists, the nearest-neighbor 
distance will be given by 



4^2 



1/3 



(48) 



The preceding analysis is based on the assumption that 
the system forms a crystal. As noted previously, the ar- 
guments based on large Nc imply that the system will 
either form a crystal or an amorphous solid. The argu- 
ment given above can serve to rule out the possibility 



of an amorphous solid. Note that the structure with the 
lowest energy is the one with the maximum atom packing 
fraction, P. Since Kepler first conjectured it in the seven- 
teenth century, it has been believed that maximal close- 
packing occurs for the HCP and FCC crystalline struc- 
tures and not for an aperiodic configuration. A recent 
proof that this is indeed the case has been constructed, 
albeit one requiring a computer to verify a large number 
of cases[T7|. 

To proceed, we need to compute the interaction en- 
ergy per baryon or, equivalently, the interaction energy 
density divided by the baryon density £\nt / p . As noted 
above, at low densities, the interaction energy between 
baryons is entirely dominated by nearest neighbor inter- 
actions. Thus, the energy per baryon from the mean-field 
energy functional Eq. ( 46 ) reduces to 



P J '^Mq 2 



J2 ffd'r d\' ^*(^1'^:(^')V'(r')0.(rO + ^*(r')0:(rOV^(r')0.(r') ^^^^ 



2 ^JJ \r-r'\ 



ip is the wave function of the quarks in one of the baryons, 
the summation is over nearest neighbors and <f>i repre- 
sents the wave function of the quarks in the nearest neigh- 
bor. We wish to minimize this, subject to constraints 
that the wave functions are orthonormal and are cen- 
tered a distance d apart. As a first step let us find the 
variational equation which minimizes £ / p subject to the 
constraint of orthonormality but without yet imposing 
the condition that the wave functions are centered a dis- 
tance d apart. A trivial application of the calculus of 
variations on £ / p with respect to ■0* yields 



2Mq j \r — r' 

.3 . <l^l{r'mr') 



U^U^+as j d^ 



(50) 



where e is a Lagrange multiplier enforcing the normaliza- 
tion of ip, and the A^ are Lagrangian multipliers enforcing 
the orthogonality of the wave function of the baryon, ip, 
and its nearest neighbors, 0;. 



Note that the terms on the left-hand side of Eq. ( 50 ) 



would be zero for a wave function that minimizes the 
single-baryon problem as in Eq. (Isl. Thus it is natural 



to treat the right-hand side of Eq. (50) as a perturba- 



tion. Treating this perturbatively has the virtue of keep- 
ing the wave functions well localized a distance d apart. 
There are two effects in this perturbation: one arising 
from the Pauli principle through the orthogonality con- 



straint encoded by the Lagrangian multipliers, and the 
other from the explicit interaction seen in the integral 
term. It is highly plausible that the second effect is para- 
metrically smaller than the first by a factor which scales 
as {asMqd)~^: the integral depends on |r — r'|~^ for wave 
functions which peak a distance d apart and which have 
a natural distance scale ~ (a^Mg)^^ . Accordingly, we 
will drop the term at this stage and subsequently verify 
that the effect is indeed small for large d {i.e., small p) 

The Lagrange multipliers A^ have dimensions of mass. 
It is convenient to rewrite them as 



A; = eoXi 



(51) 



where eo = — ■^^^js^^S.v^ ^ .16277, a'^Mg is the eigenvalue 
of the unperturbed single-baryon problem of Eq. (|8| and 
Xi is dimensionless. Since the overlaps between the wave 
functions are expected to be exponentially small, the 
are useful as expansion parameters. Thus, we can write 



i, = ij° + Y^X,^]+OiX^) 



,0 



-A,0i+O(A2) 
Y,X,e]+OiX') 



(52) 



where the sum over j indicates nearest neighbors and 
the shifts in the 4>j are only those induced by its in- 
teraction with V'- Note that at lowest order, the wave 
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functions are simply those of noninteracting baryons so 
that (f>^{r) = ip^^r — dfii) where fii is a unit vector in the 
direction between the reference particle and its nearest 
neighbor and d is the distance between nea rest neighbors 
and is well approximated by Eq. (30). It is worth 



where the wave functions are represented in Dirac nota- 
tion and p is the momentum operator. The derivation of 
this equation is straightforward but somewhat involved, 
and can be found in Appendix A. 



noting that one can only use this expansion consistently 
up to first order in the A^; higher order term will turn 
out to be parametrically smaller than the effects of next- 
to-nearest neighbor interactions which we are neglecting 
here. Inserting the expansion of Eqs. (52) into Eq. (50), 



The numerators in Eq. ( 53 ) are the squares of the rele- 



neglecting effects which are parametrically suppressed by 



vant overlap integrals Aj ^TV'^l'/'j) = / '^'^^ ^p'^{'r) ?/;°(r— 
dhj). Thus, Aj is a measure of orthogonality of two wave 
functions and and is exponentially small. Since the 
original wave functions are spherically symmetric, Aj de- 
pends only on d and not on hj — the direction connecting 



factors of (5aM,d)-i and solving for the interaction en- ^^^^ "^'.^'^e^* neighbors. Thus .4, = A for all j and all 



ergy to first order in Aj, one obtains 



terms in Eq. (53) are equal. Since A does not depend on 



mm 



2{r\ 



(53) 



the direction, without loss of generality we can take the 
direction to be along the z axis and exploit cylindrical 
symmetry: 



oo roQ 



A^27r I I dzdrr^"{y/r^ + z^)^"{y/r^ + {d~zf). 

'-00 ^0 



(54) 



r 



Moreover, it is easy to see that on the scale of d the in- 
tegrand for the r is sharply peaked and can be evaluated 
using the standard logic underlying the steepest descents 
approximation; the relative error in making this approx- 
imation is of order (figMgd)"^. The z integral can then 
be evaluated straightforwardly yielding. 



V2 



A{p)=i~p ^'^cxp(2(-)2/3. 



with p = 



and 7 « 0.044984 ; 



we have used Eq. ( 48 ) to eliminate d in favor of p. 
The numerators in Eq. (53) are given by \A\'^. 



(55) 



The 



easily evaluated: 



denominators in Eq. (53), B = 2{ips\—:2^ l^Js) , are 



B = 2 



d^p KpIiPs 



with (pIiJjs) = / d^x e-'P-^Tpsix) 



yielding B « 76.914((52Mq)-i 



(56) 



Inserting A and B into Eq. ( 53 ) and using the fact that 



all four nearest neighbor interactions are equal yields an 
interaction energy per baryon of 



SE: 



^int 
P 



.00042858 iVc Mq a/ p^.aseg 



exp(-2.0332p 



nl/3^ 



with p i 



(57) 



Equation (57) gives the interaction energy per baryon for this simplified problem. In its derivation we neglected 



the explicit potential energy contribution to the perturbation in Eq. (50) and only included the effect of the Pauli 
principle. Let us now show a posteriori that this is justified. The energy associated with this potential energy, at first 
order in the perturbation expansion, is given by 



V^.. = -m^a. I d\ d^r' ^"(^^^^'^ - Mmr^ir^ 



(58) 



r 



where n is an arbitrary unit vector and we have used the 
fact that the wave functions are real. Note that up to 
a sign, the form of this is precisely that of an ordinary 



electrostatics problem for the repulsive Coulomb energy 
of a charge distribution proportional to 'ipo{'r)'ipo{r — fid); 
the charge distribution is cylindrically symmetric and has 
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a characteristic length of d and a characteristic width of 
order (SsMq)"^. Note that the total "charge" in this 
electrostatic problem is A. It is straightforward to show 
that for asMqd > 1 that Vi^t ^ \ \og{asMQd)d-^ (up 
to overall constants). Comparing this to the expression 
for SE, one sees that 



SE 



logjasMgd) 

aU'Igd 



(59) 



which goes to zero at large d. 

It is not immediately clear whether the phase described 
in this section is absolutely stable. At the large interpar- 
ticle distances studied here, the repulsive interactions in- 
duced by the Pauli principle at the quark level dominate 
over the color-Coulomb attractive interactions from the 
Fock term. However, as the densities increase so does 
the effect of the color-Coulomb interaction. When the 
densities reach the order of (a^Mg)'^ , both the repul- 
sive effects induced by the Pauli principle at the quark 
level and the attractive Fock interactions contribute to 
the interaction energy per baryon an amount which is 
parametrically of order a^Mg. Without detailed calcu- 
lations it is impossible to know which effect is larger, i.e. 
whether the net interaction energy is positive or nega- 
tive. Note that such calculations are highly nontrivial in 
this regime since the system need not be dominated by 
nearest neighbor interactions. If it happens that interac- 
tion energy everywhere in this regime is repulsive, then 
the phase we have computed is presumably absolutely 
stable. If, however, the interaction energy is attractive 
somewhere in this regime, then the phase we have com- 
puted is not absolutely stable; the system can lower en- 
ergy by collapsing to these higher densities. However, 
even if it turns out that the low density phase is not ab- 
solutely stable in a global sense, the phase is metastable; 
all local changes in the system will necessarily raise the 
energy per baryon. 



B. The full problem 

In the last subsection we found the energy density for 
a system in which all of the quarks were constrained to 
be in the same spin-isospin state. In this section we re- 
lax this constraint to find the energy density of a sys- 
tem with fixed (low) total baryon density. Again, we 
will assume the system is in the form of a single Slater 
determinant. To create localized color singlet clusters 
(baryons) within a Slater determinant we still need to 
consider products over color of single-particle states with 
a single spin-space- flavor state and then take a product of 
these. The difference in this problem is that space- flavor 
combinations in these clusters need not be the same for 
all clusters. Thus our state is of the form 



I*)- n n 



dVV'j(r)V'sf,;a(^ 



where sij represents the spin-flavor state of the quarks 
in baryon j and a is color. 

Now let us consider the interaction between two 
baryons. At low densities the dominant interaction is 
the Pauli principle at the quark level. Note that if the 
spin-flavor state of the quarks in the two baryons are or- 
thogonal to each other the Pauli principle does not apply 
and there is no Pauli repulsion between them. If they 
are not orthogonal, there is repulsion due to the Pauli 
principle. Thus we can decrease the energy per baryon 
relative to the simplified problem of the previous subsec- 
tion by simply constructing 2Nf copies of the state from 
the simplified problem, each associated with an orthogo- 
nal spin-flavor state for the quarks: 

I*) = n n n n / \^^) ■ 

j=l a=l s = l /=1 ^ 

(61) 

Note that by construction this state is both a spin and 
flavor singlet since for each spatial wave function and 
each color one has a fully antisymmetric state in flavor 
and spin separately. The energy density is trivial to com- 
pute. As noted above there is no Pauli repulsion between 
baryons with orthogonal spin-flavor states of the quarks. 
It is also trivial to see that the color-Coulomb Fock in- 
teraction vanishes between baryons with the same spatial 
wave functions but with orthogonal spin-flavor states of 
the quarks. Thus, the total interaction energy is simply 
2Nf times the interaction energy of the simplified prob- 
lem for a simplified problem with a density of p/{2Nf). 
The interaction energy per baryon is then given by 



.00042858 Nc Mq d,^ p^-^^^^ exp ( -2.0332p 



with p : 



P 



2Nf{asMQf 



(62) 



(60) 



Note that at low densities this phase is much smaller 
than for the simplified problem. As with the simplified 
problem we do not know whether the low density phase 
described in this section is absolutely stable or merely 
metastable. However, we do know that it is at least 
met ast able. 

In concluding that this phase is at least metastable 
an issue arises which does not come up in the simplified 
problem: the role of the spin-flavor degree of freedom. 
Note that the ansatz of Eq. (61 1 has baryons which are 
not well separated. Of course, the baryons which are 
not well separated are also not interacting since they are 
made of quarks which are orthogonal in spin-flavor. How- 
ever, to conclude that the system is metastable, we need 
to know that all local variations in the state consistent 
with the Slater determinant form (which by assumption 
we take to yield the energy at leading order in the 1/Nc 
expansion) which keep the average baryon density fixed 
act to raise the energy. Clearly, in verifying this we do 
not need to consider interactions between well-separated 



16 



baryons. These act as they do in the simphfied problem. 
However, one does need to consider what happens if any 
of the 2Nf baryons with distinct spin-flavor states which 
are locahzed in the same spatial state are varied so that 
they are no longer completely orthogonal in spin-flavor 
space. Of course, if the system docs this it must simulta- 
neously alter its spatial configuration to respect the Pauli 
principle at the quark level. Clearly, the effect due to the 
Pauli principle is entirely repulsive. When the system 
makes local variations at lowest order, this Pauli effect 
is the only effect contributing. Larger scale variations 
might in principle yield attractive effects which offset the 
Pauli repulsion and lead to a net lower energy per par- 
ticle than the phase considered here. Thus, we do not 
know the system is globally stable. However, since all 
small variations only involve the repulsive interactions, 
the system is at least metastable. 



C. Subleading contributions and saturation 

So far in this section, we have shown that there ex- 
ists a low-density phase of QCD in the combined large 
Nc and \/Nc limits for which the interactions are repul- 
sive and thus a low-density phase exists which is at least 
metastable. Such a phase clearly does not lead to sat- 
urating matter in the manner of QCD at Nc = 3 and 
physical quark masses in which the matter is self-bound 
and stable with nonzero density and zero external pres- 
sure. We do not know at this stage whether the phase we 
have studied is absolutely stable in the combined limit or 
whether there are attractive interactions at shorter dis- 
tances that might give rise to a saturating phase. 

Remarkably, however, it is possible to conclude that 
as the combined large Nc and heavy quark limit is ap- 
proached, the system must have a saturating phase pro- 
vided that the lightest glueball in the spectrum at large 
Nc is scalar. 

At first blush, it may seem that, by construction, sub- 
leading effects in 1/Nc must have a small effect on the in- 
teraction energy and thus as the limit is approached they 
are extremely unlikely to affect qualitative issues such as 
saturation. This is not correct, however, since the prob- 
lem involves multiple limits — the low density limit and 
the large Nc limit and the heavy quark limit — and these 
need not commute. Note that in the analysis done above 
we have first gone to the combined heavy quark and large 
Nc limits at fixed p/ {Mq ds)^ and then considered the 
low density behavior. The interaction energy per baryon 
is then given parametrically by Nc Mq dg"^ while its den- 
sity dependence is of the form / [p/{Mq dsY^ where / 
goes to zero with p. Suppose that there is an additional 
contribution to the interaction energy baryon which is 
parametrically of the form 



5E, 



■gb 



AqcdA^°«.' /gb (M^cd) 



(63) 



with /gb going to zero with p. In the limit in which we 
worked, with heavy quark and large Nc limits at fixed 



p/{Mq ds)^ taken prior to the low density limit, this 
contribution is clearly subleading. However, if the low- 
density limit is taken first SEg^ becomes the dominant 
term. 

If SEgh is attractive at the lowest densities, the system 
at very low density will be unstable, ultimately collapsing 
to a density in which its attraction is counteracted by 
repulsive interactions 6E (at leading order in Nc yielding 
saturated matter. Note that as the large Nc and heavy 
quark limits are approached, the density at which the 
pressure of the repulsive leading-order interaction cancels 
out the attraction in SEgy, is pushed towards zero as is the 
interaction energy. Thus, there will be a very low-density 
regime in which there is saturating matter, a low-density 
repulsive regime where our previous calculations apply, 
and a higher-density regime where we do not presently 
know the equation of state of cold matter. 

The preceding scenario depends on the existence of an 
attractive interaction energy of the form of Eq. ( 63 1 at 
the lowest of densities. An equivalent way to state this 
is that since the mass of the baryons are heavy in our 
limit, baryon kinetic energies are suppressed so that the 
scenario depends on the inter-baryon potential between 
baryons having as its longest range contribution some- 
thing attractive with a range of order Aqcd- In this 
combined limit the longest-range interaction is due to 
glueball exchange — meson exchange has a range of char- 
acteristic order of Mq^ and Pauli effects at the quark 
level have a range {dsMQ)~^. Note that a one-glueball 
exchange potential has a strength which is parametrically 
of order Nc in the 1/A^c expansion (which is parametri- 
cally of relative order 1/Nc down compared to the lead- 
ing order contributions). Moreover in the heavy quark 
limit, the baryon is small and thus gluonic couplings are 
expected to be small. One expects that the dominant 
coupling of a glueball to the quarks in the baryon will be 
via two gluons so that the glueball-baryon coupling will 
scale as ds and the one-glueball exchange potential will 
scale as d^. This characteristic Nc and dg dependence is 
seen in Eq. (63). Now suppose that the lightest glueball 
at large Nc is a positive parity scalar — as we expect — 
the longest-range interaction is attractive and the system 
saturates with a saturation density parametrically small 
compared to (Mga^)'^ (the natural length scale of baryon 
interaction) . Thus we expect that the saturation density 
will go and one which goes to zero as the combined large 
Nc and heavy quark limit is approached. 

One feature might simplify the computation of the sat- 
uration density due to subleading effects: the fact that 
the saturation density is pushed to zero as the combined 
limit is approached. This in turn means that only the 
longest-range interaction due to glueball exchange — i.e., 
the exchange of the lightest glueball — is relevant. 

However, the problem of computing the saturation 
density and energy density has some non-trivial features. 
Note that in the purely repulsive case studied in the pre- 
vious section, the system was in a close packed crystal in 
order to minimize the repulsive interactions. Once one 
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includes the attractive non-leading but long-range effects 
due to glueball exchange it is not obvious a priori which 
crystal structure is optimal. Moreover, the spin-flavor 
effects now come into play in a more consequential way, 
Note that at leading order in 1/Nc, baryons composed of 
quarks of orthogonal spin-flavor states did not interact. 
However, they do feel one-glueball exchange. There is 
an additional complication: in order to make a detailed 
calculation we need to start with a reliable energy func- 
tional. We believe that we have a reliable energy func- 
tional for inter-baryon distances of order MqUs but that 
is numerically small when working in the combined limit. 
We also believe that a one-glueball exchange potential is 
valid when working at very long distances. However, to 
find the saturation density we are necessarily in a regime 
where the two effects yield equal and opposite pressures. 
We could presumably attempt to compute this by taking 
as our energy functional the sum of the one from Pauli 
repulsion and the one-glueball exchange. However, from 
a theoretical perspective, it is not completely clear that 
such a prescription is valid for this regime. 

Given these complexities we will defer any quantitative 
study of this problem to the future. At this stage, we 
merely note the qualitative result that the system will 
saturate. 



IV. DISCUSSION AND OUTLOOK 

In this work we have explored "nuclear physics" in an 
an artificial world in which the number of colors is large, 
as is the ratio of the quark masses to Aqcd- We have 
computed properties of both baryons and cold baryon 
matter in the combined limits. While we have made 
significant progress in learning about this world, some 
important open problems remain. 

The baryon's mass and form factor were computed nu- 
merically to high accuracy in the combined limit using 
the mean-field formulation of WittenfJ. The problem at 
leading order is essentially solved. However, there is an 
important open problem concerning corrections to the 
leading order result. 



HE the 



As discussed in Subsec. 
mean-field calculation correctly captures the leading or- 
der of the 1 /Nc expansion (with the heavy quark limit im- 
plicitly taken at the outset). However, it was also shown 
there that sub-leading terms in the 1/Nc expansion does 
not correspond to perturbation theory in the residual in- 
teraction beyond the one-body mean-field Hamiltonian. 
This raises an important challenge; namely, to formulate 
a tractable scheme for the systematic inclusion of sub- 
leading terms in the 1/Nc expansion. Of course, even if 
one can compute such corrections systematically there is 
still the problem of how to reliably compute corrections 
to the heavy quark limit. 

Regarding baryon matter, we have derived an expres- 
sion for the energy density at leading order in the com- 
bined expansion for a phase of matter at fixed density 
which is taken to be on the scale of (Mqas)^ but numer- 



ically small on this scale. This phase is repulsive and has 
an exponentially small energy density. We have shown 
that it is at least metastable. We have also shown that 
subleading effects lead to saturating matter in which the 
saturation density is parametrically small on the scale 
of (Mqcis)^ in which the saturation density and energy 
density both tend to zero as the combined limit is ap- 
proached. However there are a large number of open 
problems concerning baryon matter. 

At a theoretical level, the most basic issue concerns 
the validity of the mean- field treatment. In this work, 
we assumed that a mean-field theory treatment based on 
a single Slater determinant was valid in the combined 
limits. This assumption was implicitly used in ref. [7| in 
its treatment of problems with more than one baryon. 
As we remarked earlier, while this assumption is plausi- 
ble, it has not been demonstrated with anything like the 
degree of rigor for the single baryon problem. Thus, the 
most critical outstanding theoretical problem is to deter- 
mine whether or not this assumption is valid. A related 
problem concerns corrections to the mean-field result. As 
in the single baryon case, we have no tractable scheme 
to include 1/Nc corrections and corrections to the heavy 



quark limit. As noted in Subsec. HI C a viable scheme to 



reliably include the longest distance effects (even though 
they are subleading) is essential if one wishes to com- 
pute the properties of saturating nuclear matter which 
arises from the interplay of subleading but long distance 
attractive effects of glueball exchange with the leading 
order but shorter-range repulsion arising from the Pauli 
principle. 

Apart from these theoretical issues, it is worth noting 
that the phenomenology of baryonic matter in this 
artificial world is not well understood except at low 
density, even in the combined limit. Assuming that the 
mean-field treatment can be justified, in principle one 
could study the problem of baryonic matter numeri- 
cally and thereby access higher densities. Perhaps the 
most important question that needs to be addressed 
is whether baryonic matter saturates with densities 
on a scale of {Mqas)^- If it does, such a phase will 
prove to be the true ground state of baryonic matter 
in this regime rather than the phase arising from the 
interplay of subleading, long distance attractive effects 
of glueball-exchange with the leading order but shorter- 
range repulsion. A final important question is whether 
baryonic matter undergoes any phase transitions as the 
density is increased. 

This work is supported by the U. S. Department of 
Energy through grant number DEFG02-93ER-40762. 



Appendix A: The interaction energy 



This appendix contains the derivation of Eq. (53). 



We begin by inserting the first-order expansions of -0*, 
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ip, and e in A; Eq. (50) becomes 
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Equating like powers of \j yields: 
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The zeroth-order equation is simply the Schrodinger 
equation for the isolated baryon. It is useful to rewrite 
the first order equation in Dirac notation, 
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^)|^j)-e]|^°)=e|0°) 



where Hq is the one-body mean-field Hamiltonian for the 
isolated baryon. Acting on both sides with the bra 
yields 

-e] = e'^A withyt= (V^°|0°). 

where by construction, the value of A is independent of 
j since all the nearest neighbors are equivalent. 

The interaction energy at leading order depends on 

The next step is 



Aj-e^. We have already computed e 



to compute \j . To proceed, we insert into Eq. ( A2 1 , 
yielding 



2M. 



+ Vo-eA \^]) = eo (|(/),) - Am) . (A3) 



inverse operator is well defined on the space of square 
integrable functions since is a positive operator while 
eo is negative and the operator is thus finite acting on 
any wave function. Note, moreover, that Eq. (A4| does 



indeed yield a wave function peaked around r — dfij with 
a characteristic size of order {asMQ)~^\ this justifies the 
approximation. 

Recall that Aj enforced the orthogonality of -0 and 4)j : 

((V/H-A,(^]|)(|0°) + A,|</-]))=O 
^ + A,((V/|0,") + (0"|(/.])) =0 
A + 2Xj{i,^\(l)°) 



where the last equality follows by symmetry. We expect 
that the effect of the Pauli principle on each baryon in a 
pair of nearest neighbors to be identical. Using this last 
form and our expression for j-j/'j) yields 



On physical grounds, we expect to be peaked 

around r — dfij with a characteristic size of order 
((5sMq)~^ since that is the type of wave function which 
will efficiently yield orthogonality between ip and 4>j . Let 
us start by doing the calculation assuming this to be true 
and then verify a posteriori that the assumption was cor- 
rect. With this assumption it is easy to see that the 
contribution to the left-hand side of the Vq term is char- 
acteristically of order ag/d while the eo term is of order 
alMq. Thus, for low densities (large Mq d ^ 1) one 
can drop the contribution of Vq so that 
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where in the last form if)'^ in the denominator has 
been replaced with ^p'^. This is legitimate since the 
operator( ^^^^ — e°)~^ is translationally invariant and the 
potential term is not present; thus two expressions are 
equivalent. 

The total shift in energy per baryon is thus given by: 



2m 



(A5) 



where we have also dropped the term proportional to 
A since it is exponentially suppressed. Note that the 



which is the form of Eq. ( 53 I . 
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